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We introduce a natural language interface for building stochastic % calculus models of biological 
systems. In this language, complex constructs describing biochemical events are built from basic 
primitives of association, dissociation and transformation. This language thus allows us to model 
biochemical systems modularly by describing their dynamics in a narrative-style language, while 
making amendments, refinements and extensions on the models easy. We demonstrate the language 
on a model of Fey receptor phosphorylation during phagocytosis. We provide a tool implementation 
of the translation into a stochastic % calculus language, Microsoft Research's SPiM. 1 2 



1 Introduction 

Modelling of biological systems by mathematical and computational techniques is becoming increas- 
ingly widespread in research on biological systems. In recent years, pioneered by Regev and Shapiro's 
seminal work lTT6l[T7l . there has been a considerable amount of research on applying computer science 
technologies to modelling biological systems. Along these lines, various languages with stochastic sim- 
ulation capabilities based on, for example, process algebra, term rewriting (see, e.g, 13 H) and Petri 
nets (see, e.g., |fl9l l9ll) have been proposed. However, expressing biological knowledge in specialised 
modelling languages often requires a simultaneous understanding of the biological system and expert 
knowledge of the modelling language. Isolating and communicating the biological knowledge to build 
models for simulation and analysis is a challenging task both for wet-lab biologists and modellers. 

Writing programs in simulation languages requires specialised training, and it is difficult even for the 
experts when complex interactions between biochemical species in biological systems are considered: 
the representation of different states of a biochemical species with respect to all its interaction capabilities 
results in an exponential blow up in the number of states. For example, when a protein with n different 
interaction sites is being modelled, this results in 2" states, which needs to be represented in the model. 
Enumerating all these states by hand, also without inserting typos is a difficult task. 

*This work has been done during Kahramanogullari's appointment at the Department of Computing, Imperial College and 
Centre for Integrative Systems Biology at Imperial College. Kahramanogullan acknowledges support of the UK Biotechnol- 
ogy and Biological Sciences Research Council through the Centre for Integrative Systems Biology at Imperial College (grant 
BB/C5 19670/1). 
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To this end, we introduce an intuitive front-end interface language for building process algebra mod- 
els of biological systems: process algebras are languages that have originally been designed to formally 
describe complex reactive computer systems. Due to the resemblance between these computer systems 
and biological systems, process algebra have been recently used to model biological systems. An impor- 
tant feature of the process algebra languages is the possibility to describe the components of a system 
separately and observe the emergent behaviour from the interactions of the components (see, e.g., HII3)- 

Our focus here is on % calculus lPT2l . which is a broadly studied process algebra, also because of 
its compactness, generality, and flexibility. Because biological systems are typically highly complex 
and massively parallel, the % calculus is well suited to describe their dynamics. In particular, it allows 
the components of a biological system to be modelled independently, rather than modelling individual 
reactions. This allows large models to be constructed by composition of simple components. % calculus 
also enjoys an expressive power in the setting of biological models that exceeds, e.g., Petri nets 0. 

In the following, we present a language that consists of basic primitives of association, dissociation 
and transformation. We impose certain consistency constraints on these primitive expressions, which 
are required for the models that describe the dynamics of biochemical processes. We give a translation 
algorithm into stochastic % calculus. Based on this, we present the implementation of a tool for automated 
translation of models into Microsoft Research's stochastic simulation language SPiM lfl4l[T3l . which can 
be used to run stochastic simulations on % calculus models. We demonstrate the language on a model of 
Fey receptor phosphorylation during phagocytosis. The implementation of the translation tool as well as 
further information is available for download at our website 

2 Species, Sites, Sentences and Models 

We follow the abstraction of biochemical species as stateful entities with connectivity interfaces to other 
species: a species can have a number of sites through which it interacts with other species, and changes 
its state as a result of the interactions (HOD]]. m Section |4j we use this idea to design a natural language- 
like syntax for building models. The models written in this language are then automatically translated 
into a SPiM program by a translation algorithm. This is done by mapping the sentences of the lan- 
guage into events constructed from the basic primitives, which are then compiled into executable process 
expressions in the SPiM language. 

There is a countable set of species A,B,C, Each species has a number of sites a,b,c,... with 

which it can bind to other species or unbind from other species when they are already bound. We write 
sentences that describe the 'behaviour' of each species with respect to their sites. There are three kinds 
of sentences, i.e., associations, dissociations, and transformations. We define the sentences as 

(type, (A, a), (B,b), Pos, Neg, r) 

where type € {association, dissociation, transformation] is the type of the sentence. The pairs (A, a) 
and (B,b) are called the body of the sentence. The sets Pos and Neg are called the conditions of the 
sentences. (A, a) and (B,b) are pairs of species and sites, and Pos and Neg are sets of such pairs of 
species and sites. If the sentence is an association, it describes the event where the site a on species A 
associates to the site b on species B if the sites on species in Pos are already bound and those in Neg are 
already unbound. If it is a dissociation sentence, it describes the dissociation of the site a on species A 
from the site b on species B. A transformation sentence describes the event of species A transforming into 
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species B, where B can be empty. In this case, this describes the decay of species A. In transformation 
sentences, sites a and b must be empty, since transformations are site independent, r G M + denotes the 
rate of the event that the sentence describes. Then a model ^# is a set of such sentences. In Section 
|4| we give a representation of these sentences in natural-language. For example, a sentence of the form 
(association, (A, a), (B,b), {(A,c)}, { }, 1.0) is given with the following English sentence. 

site a on A associates site b on B with rate 1.0 if site c on A is bound 

We denote with species(^#) all the species occurring in the body of the sentences of jM. The func- 
tion sites(^#,A) denotes the sites of the species A that occur in the body of all the sentences of j$ . 
sites^os, A) denotes the sites of the species A in Pos (similarly for Neg). For any set £/, p(stf) denotes 
the powerset of srf ' . 

2.1 Conditions on Sentences 

Given a model jM, we impose several conditions on its sentences. 

1. Sentences contain relevant species. The species in the condition of each sentence must be a 
subset of those in the body of the sentence. 

2. Conditions of the sentences are consistent. For every sentence of the form 
(type, (A, a), (B,b), Pos, Neg, r), we have that PosHNeg = 0. 

3. All the sites in the conditions are declared in the model. For every sentence of the form 
(type, (A, a), (B,b), Pos, Neg, r), we have that sites(Pos,A) C sites(^#,A), 

s\tes(Neg,A) C sites(^#,A), s\tes(Pos,B) C sites(^#,B) and s\tes(Neg,B) C sites(^f 

4. Association sentences associate unbound species. For every association sentence 
(association, (A, a), (B,b), Pos, Neg, r), we have that (A, a), (B,b) G Neg. 

5. Dissociation sentences dissociate bound species. For every dissociation sentence 
(dissociation, (A, a), (B,b), Pos, Neg, r), we have that (A, a), (B,b) G Pos. 

6. Transformation sentences are unbound at all sites. For every transformation sentence 

( transformation^, B, Pos, Neg, r), we have that Pos = and Neg = {(A,x) \ x G sites(^#,A)}. 

When these conditions hold, we can map the sentences of a model to another representation where the 
role of the conditions become more explicit. In the following, for a model we describe the states of its 
species as subsets of its sites, where bound sites are included in the set describing the state. For example, 
for a species A with binding sites sites(^#,A) = {a,b}, the set p(s\tes(^ ,A)) = {{},{a},{b},{a,b}} 
is the set of all its states. Then {a} is the state where site a on A is bound and site b on A is unbound. We 
map each sentence (type, (A, a), (B,b),Pos, Neg, r) to a sentence of the form 

(type, (A, a), (B,b), states(A), states(B), r) 

where states(A) and states(B) are obtained as follows. 

states(A) = {,¥ G ^(sites(^,A)) | ((A,x) G Pos =>- x G y) A (x G y =>- (A,x) <£Neg)} 

This representation allows us to impose another condition on the sentences: 
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7. There are no overlapping conditions in the sentences. For any two sentences of a model ^# 
of the form (type 1; (A, a), (B,b), Pos\, Negy, r) and (type 2 , {A, a), (B,b), P0S2, Neg 2 , r) where 
typej = type 2 , we obtain states(A)i and states(fi)i, for the first and states(A)2 and states(5)2, 
for the second sentence. Then we have that 

- if states(A)i = states(A)2 then it must be that states(B)i n states(B)2 = 0; 

- if states(fi)i = states(B) 2 then it must be that states(A)i n states(A) 2 = 0; 

- if states(A)i 7^ states(A)2 and states(Z?)i 7^ states(B)2 then it must be that 
states(A)i n states(A) 2 = and states(fi)i n states(B) 2 = 0. 

Example 1 Consider the models 



Jt\ = { (association, (A,a), (B,b), {(B,f)}, {(C,c), (A, a), (B,f)}, 1.0), 
(dissociation, (A, a), (B,b), {(B,b)}, {}, 1.0), 
(transformation, A, B, {}, {}, 1.0), 
(association, (D,d), (E,e), {}, {(D,d), (E,e)}, 2.0), 
(association, (D,d), (E,e), {}, {(D,d), (E,e)}, 4.0)} 



This model does not fulfill any of the conditions above: in the first sentence, (I.) C ^ {A,B}; (2.) 
(BJ) £ Pos and (B,f) G Neg; (3.) f £ {b}; (4.) (B,b) £ {(C,c), (A,a),(B,f)}. In the second sentence, 
(5.) (A, a) {(B,b)}. In the third sentence, (6.) {} 7^ {(A, a)}. (7.) In the fourth and fifth sentences, 
states(D)i = {{}} = states(D)2 and states(£)i = {{}} = states(£")2. 

Example 2 The model ^2 fulfills all the conditions above. 



3 Translation into Stochastic n calculus 

We use the representation of the states of species as sets of their sites to map models to stochastic % 
calculus specifications. For this purpose, we first map each model to a compile map. Let us first recall 
some of the definitions of stochastic % calculus, implemented in S Pi M, as they can be found in |[T3ll . 

3.1 Stochastic % calculus 

Definition 3 Syntax of stochastic K calculus. Below fn(P) denotes the set of names that are free in P. 



^#2 = { (association, (A,a\), (B,b), {}, {(A,a\), (B,b)}, 1.0), 
(association, (A, aj), (C,c), {}, {(A, 02), (C,c)}, 1.0), 
(dissociation, (A,ai), (B,b), {(A,a\), (A, 02), (B,b)}, {}}, 2.0). 
(dissociation, (A,a\), (B,b), {(A,ai), (B,b)}, {(A,a 2 )}}, 4.0) } 
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Expressions above are considered equivalent up to the least congruence relation given by the equivalence 
relation = defined as follows. 

P|0 = P newxO = 

P\Q = Q\P newi iiew}'f = new^newif 

P\{Q\R) = (P\Q)\R newx(P|£) = P \ new* Q ifx $jh(P) 

X(n) = P {n/m} ifX(m)=P 

3.2 Compile Maps 

We map models into compile maps, denoted with % '. A compile map is a set of expressions that we 
call process descriptions for each species A G species(^#). For a model the process description 
of species A G species(^), denoted with P(A), is the pair (A, actions(A)). Here, actions(A) is the set 
collecting actions(A,^) for every y £ p(sites(^#,A)). 

actions(A,^) = (^,assoc(A,^),dissoc(A^),transform(A,^)) 

We define assoc(A, I 5^) as the set of assoc(A,,5^,a) for every a £ sites(^#,A). 

assoc(A,^,a) = (a,assocPartners(A,^,a)) 

where assocPartners(A, I 5^,a) is the set 

{(B,b,states(B),r) | 

( {association, (A, a), (B,b),Pos, Neg, r) £ J% A 5? £ states(A) ) 

V ((association, (B,b), (A, a), Pos, Neg, r) £ Jt A 5? G states(A))}. 

We define d\ssoc(A,y), similarly, as the set of dissoc(A, t 5^,a) for every a £ sites(^#,A). 

dissoc(A,^,fl:) = (ajdissocPartners^,^,^)) 

where dissocPartners(A,^,a) is the set 

{(B,b,states(B),r) \ 

((dissociation, (A, a), (B,b), Pos, Neg, r) £ ^# A y £ states(A) ) 

V ((dissociation, (B,b), (A, a), Pos, Neg, r) £ Jt A £ states(A) ) } . 

If 5? = 0, the set transform(A, o5^) is defined as 

{ (B, r) | ( ( transformation, A, B, Pos, Neg, r ) £ } . 

Otherwise, it is 0. 

Example 4 Consider the model ^M 2 m Example^ We have that the compile map ¥? 2 for this model is as 
follows. 

{ (A, { ({}, {( ai ,{(B,b,{{}}, 1.0)}), (a 2 ,{(C,c,{{}},\.0)}, {}, {}), 

{{ ai },{(a 2 ,{(C,c,{{}},l.0)},{(B,b,{{b}},4.0)},{}), 

({a 2 },{( ai ,{(B,b,{{}},l.O)}, {},{}), 

{{ ai ,a 2 },{},{(B,b,{{b}},2.0)}, {})}), 
(B,{ ({},{(b,{(A, ai ,{{},{a 2 }},\.0)}, {},{}), 

{ {b}, {}, {(b,{(A, ai ,{{ ai }}, 4.0), (A, ai ,{{ ai ,a 2 }}, 2.0)}, {})}), 
(C,{ ({},{(a l ,{(A,a 2 ,{{},{a l }},l.O) },{},{}>, 

(H, O.O. {}>}» 
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3.3 From Compile Maps to Stochastic n calculus 

We construct a % calculus specification from the compile map 'if of a. model For each species 
A G species(^), we map the process description P(A) to a process specification in stochastic % calculus. 
Let 

P(A) = (A,{actions(A,^i),...,actions(A,^„)}) 

where p(sites(^#,A)) = {S?\ , . . . , J^ n }, that is, the powerset of set of sites of A. Thus, there are n process 
specifications for the species A, some of which may be empty. Each process specification for each state 
y of A is of the following syntactic form. 

process declaration "= (" local channel declarations 

association specifications 
"+" dissociation specifications 
"+" transformation specifications ")" 

The idea here is that each set of sites of a species A denotes the state where the sites in the set are 
bound. Thus the powerset of the set of sites of a species denotes the set of all its states. Now, let us obtain 
the process expression for each state ^ with respect to actions(A, = 5^) where 1 < i < n. Let us consider 

= {ai, . . . , at} of A with 

actions(A,^) = (^,assoc(A,^),dissoc(A,^),transform(A,^)) . 
Process declaration 

The expression for process declaration is a process name with its list of parameters. It is delivered by 
the dissociation sentences in ^# and S?i = {a\ , . . . , a k }. For every aj G J^-, consider the set 

&(A,aj) = 

{(aj, (r/2)) | (dissociation, (A, aj),(B,b), Pos, Neg, r) G U 
{(aj, (r/2)) | (dissociation, (B,b), (A, aj), Pos, Neg, r) G U 
{(aj, 1.0) | (dissociation, (B,b),(A,aj), Pos, Neg, r) ^ ^# A 

(dissociation, (A,aj) , (B, b) , Pos, Neg, r) ^ . 

We associate each element of the set M(A,aj) a unique label s G N + and obtain M'(A,aj). Then if 
£%'(A,aj) = {(aj, n, 1), . . . , (aj, r£, £) } we write the process declaration for A at state ^ = {a\,. . . ,a k } 
as follows. 

Ai(a\l,...,a\l\, ,a k l,...,a k £ k ) 

Example 5 For the state ^2 = { a i} of species A of Example^ we have the process declaration below, 
since we have that M'(A,a\) = {(a\,2.Q, 1), (ai, 1.0,2)}. 

A 2 (a\\,a\_2) 

Local channel declarations 



These expressions are delivered by the dissociation sentences in M and the assoc(A,^). That is, for 
every 

assoc(A,^,ay) = (aj, assocPartners(A, t 5^-,a y -) ) G assoc(A,^-) , 
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and for every (B,b, states (B), r) G assocPartners(A,^-,ay) consider the set 

W(A,aj,B,b) = 

{(ay, (r/2)) | (dissociation, (A,aj),(B,b), Pos, Neg, r) G M A ay -< b} U 
{(ay, (r/2)) | (dissociation, (B,b),(A,aj), Pos, Neg, r) G A ay -< U 
{(ay, 1.0) | (dissociation, (B,b),(A,aj), Pos, Neg, r) ^ ^# A 

(dissociation, (A, ay) , (B, b) , Pos, Neg, r) g' Aaj ~< b} 

where -< denotes a lexicographic order on sites. We associate each element of the set (A,aj,B,b) a 
unique label s G N + to obtain W(A,aj,B,b). Then if 

q/ l (A,a j ,B,b) = {(a j ,r l ,\),...,(a j , n ,l)} 

then we write the channel declarations for assoc(A,^-,ay) as follows. 

new ayl@ri ... new ajlQrf 

Example 6 For the state 5^2 = { fl i} of species A of Example^ we have the channel declarations below, 
since we have that e %?'(A,a2,B,b) = {(a2, 1.0, 1)}. 

new a2l@1.0 



Association specifications 

The expression for association specifications for species A at state assoc(A, 5^\) is delivered by assoc(A, 5f\). 
For every 

(aj, assocPartners(A, t 5^-,a ; ) ) G assoc(A,^), 
and for every (B,&,states(fi), r) G assocPartners(A,^-,a ; ) consider the set 

m(A,aj,B,b) = 

{(\ajb, r) \((B,b),states(B),r) G assocPartners(A,^-,ay) Aay-<&} U 
{(?baj,r) | {(B,b), states(B), r) G assocPartners(A,^-, ay) Ab-<aj}. 

We associate each element of the set S3(A,aj,B,b) a unique label s G N + and obtain ' (A,aj,B,b). 
Association of site ay on A results in the state = t 5^-U{ay}. For each element of (\ajb,r s ,s) G 
£$'(A,aj,B,b) we write the following, composed by +. 

! ajbs(aj 1 , . . . , ay^) ; continuation 

The association channel names, such as ajbs here, are also declared as global channel declarations, 
preceding all the process declarations. The continuation is written for A in S^t as for process declarations 
above, however we write nil for the channel names for those associations of site ay on A with some 
site b' / b. Here, nil is the nil-dissociation channel with rate 0. We obtain a/1,. . . ,aj£ from the set 
& (A,aj,B,b) as in channel declarations. 

Example 7 For the state 5^2 = {ai} of species A of Example^ we have the following association 
specifications. 

\a7.c\{a.2) ;A3(ai \,a\l,a2) 
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Dissociation specifications 

The expression for dissociation specifications for species A at state assoc(A, S^i) is delivered by dissoc(A, 
For every 

(aj, dissocPartners(A, t 5^-,a ; ) ) G dissoc(A,^), 
and for every (B,b, states (B), r) G dissocPartners(A,^-,a 7 ) consider the set 

<g(A,aj,B,b) = 

{{\aj, r) | ((B,b), states(S), r) G dissocPartners(A,^,a/) A aj -< b} U 
{(?b, r) | ((B,b), states(fl), r) G dissocPartners(A,^-,a i ) Ab^aj}. 

We associate each element of the set W(A,aj,B,b) a unique label s G N + and obtain W(A,aj,B,b). 
Dissociation of aj on A results in the stated = y t \{aj}. For each ( \aj, r s , s) £W(A,aj,B,b) we write 
the following, composed by "+". 

! ajS ; continuation + ?ajs ; continuation 

The continuation is written for A in as for process declarations above. 

Example 8 For the state = {ai} of species A of Example^ we have the following dissociation 
specifications. 

!ail;Ai() + ?flil;Ai() 

Transformation specifications 

The expression for transformation specifications for species A is given only if the state S" = {}. In that 
case, for transfrom(A, {}) = {(fii , n), . . . , (B^,r^)} we write 

delay@ri;5i() + ... + delay*/* ;B k () 



4 Syntax of the Language 

The syntax of the language is defined in BNF notation, where optional elements are enclosed in braces 
as {Optional}. A model (description) consists of sentences of the following form. 
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Model 
Sentence 



Association 

Dissociation 

Phosphorylation 

Dephosphorylation 

Transformation 
Decay 
Conditions 

Condition 

Site 

Species 



Sentence y . . . Sentence m m > 1 

Association 
Dissociation 
Transformation 
Decay 

Phosphorylation 
Dephosphorylation 

::= Site on Species associates Site on Species 

{with rate Float} {if Conditions} 

::= Site on Species dissociates Site on Species 

{with rate Float} {if Conditions} 

::= Site on Species gets phosphorylated 

{with rate Float} {if Conditions} 

::= Site on Species gets dephosphorylated 

{with rate Float} {if Conditions} 
Species becomes Species {with rate Float} 

Species decays {with rate Float} 

Condition 

Condition and Conditions 

Site on Species is bound 
Site on Species is unbound 

String 

String 



In our implementation of the translation algorithm, each sentence of a model given in this syntax is 
mapped by a lexer and parser to a data structure of the form given in Section[2]in the obvious way. Phos- 
phorylation sentences are treated as association sentences where the second species is by default Phosph 
with the binding site phosph. The dephosphorylation sentences are mapped similarly to dissociation 
sentences. If not given, a default rate (1.0) is assigned to sentences. 



4.1 A Model of Fey Receptor-mediated Phagocytosis 

We demonstrate the use of the language on a model of Fey receptor (FcyR) phosphorylation during 
phagocytosis, where the binding of complexed immunoglobulins G (IgG) to FcyR triggers a signalling 
cascade that leads to actin-driven particle engulfment [8j Q~8j SI- When a small particle is coated (op- 
sonised) with IgG, the Fc regions of the IgG molecules can bind to FcyRs in the plasma membrane and 
initiate a phagocytic response: a signalling cascade then drives the remodelling of the actin cytoskele- 
ton close to the membrane. This results in cup-shaped folds of plasma membrane that extend outwards 
around the internalised particle and eventually close into a plasma membrane-derived phagosome. 

FcyR contains within its cytoplasmic tail an immunoreceptor tyrosine-based activation motif (ITAM). 
The association of FcyR with an IgG induces the phosphorylation of two tyrosine residues within the 
ITAM domain by Src-family kinases. The phosphorylated ITAM domain then recruits Syk kinase, which 
propagates the signal further to downstream effectors (see Figure [j}. In our language, we can describe 
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Figure 1 : A simple model of the phosphorylation of the ITAM domain on the Fey receptor during phago- 
cytosis. Adapted from (S|. 

the initial phases of this cascade as follows: 

site f on FcR associates site i on IgG with rate 2.0 

site y on FcR gets phosphorylated if site f on FcR is bound 

site z on FcR gets phosphorylated if site f on FcR is bound 

The first sentence above describes the binding of FcyR to IgG. The second and third sentences de- 
scribe the phosphorylation of the two tyrosine residues on ITAM (association of a PhosphO () molecule). 
This is automatically translated by our tool into the SPiM program given in Appendix A. We can then 
run stochastic simulations on the model given by these sentences. 

By using this language and our translation tool, we can build models of different size and complexity, 
and modify and extend these models with respect to the knowledge in hand on the different sites and 
interaction capabilities of the FcyR, as well as other biological systems. For example, the model above 
abstracts away from the role played by the Src kinases in the phosphorylation of the FcyR as depicted in 
Figure [T] The sentences above can be easily modified and extended to capture this aspect in the model 
as follows: here, the shaded part demonstrates the modifications with respect to the model given above. 

site f on FcR associates site i on IgG with rate 2.0 

site y on FcR gets phosphorylated if site s on FcR is bound 

site z on FcR gets phosphorylated if site s on FcR is bound 

site s on FcR associates site sr on Src if site f on FcR is bound 

site s on FcR dissociates site sr on Src 

The SPiM program resulting from automated translation of this model is given in Appendix B. It is 
important to note that, because FcR has 4 binding sites in the model above, in the SPiM code resulting 
from the translation, there are 16 species for FcR, denoting its different possible states. However, in the 
code given in Appendix A, there are 8 species for FcR denoting its states that result from its 3 binding 
sites in that model. 

5 Discussion 

We have introduced a natural language interface for building stochastic % calculus models of biological 
systems. The fc-calculus (21 [7] and the work on Beta-binders in ifTTTl have been source of inspiration 
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for this language. 

In ifTTIl . Guerriero et al. give a narrative style interface for the process algebra Beta-binders for a 
rich biological language. In our language, we build complex events such as phosphorylation and dephos- 
phorylation of sites as instances of basic primitives of association, dissociation and transformation. We 
give a functional translation algorithm for our translation into stochastic % calculus. The conditions that 
we impose on the models are automatically verified in the implementation of our tool. These conditions 
should be instrumental for 'debugging' purposes while building increasingly large models. 

The implicit semantics of our language, which is implemented in the translation algorithm into % 
calculus, can be seen as a translation of a fragment of the K calculus into % calculus. Another approach 
similar to the one in this paper is the work by Laneve et al. in [ 15], where the authors give an encoding 
of nano-K-calculus in SPiM. In comparison with our algorithm, the encoding in [15] covers a larger part 
of nano-fc by using the SPiM language as a programming language for implementing a notion of term 
rewriting, where there is an explicit function for matching. The algorithm gives the different states of a 
species in the SPiM encoding with respect to the parameters of that species as in K"-calculus. 

Topics of future work include an exploration of the expressive power of the association, dissociation 
and transformation primitives with respect to Kohn diagram representation [ 10] of biological models. 
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Kahramanogullan, Cardelli & Caron 

Appendix A 

site f on FcR associates site i on IgG with rate 2.0 

site y on FcR gets phosphorylated if site f on FcR is bound 

site z on FcR gets phosphorylated if site f on FcR is bound 

The SPiM code resulting from the automated translation of this model. 



directive sample 10.0 
directive plot FcR7() ; FcR6(); 

FcR5(); FcR4(); FcR3(); 

FcR2(); FcRlO; 

FcROO; IgGK); IgGOO; 

PhosphK); PhosphOO 

new f il@l . : chan(chan) 
new phosphy2@l . : chan(chan) 
new phosphz3@l . : chan(chan) 
new nil@0.0:chan 

let FcROO = 

( new f@1.0:chan 

!fil(f)*2.0; FcRl(f) ) 



and FcR5(f : chan, z : chan) = 

( ?phosphy2(y) ; FcR7(f,y,z) ) 

and FcR6 (y : chan , z : chan) = 
( new f@1.0:chan 

!fil(f)*2.0; FcR7(f,y,z) ) 

and FcR7(f : chan, y: chan, z : chan) = 


let IgGOO = 

( ?fil(i); IgGl(i) ) 

and IgGl(i:chan) = 




and FcRl(f : chan) = 

( do ?phosphy2(y) ; FcR4(f,y) 
or ?phosphz3(z) ; FcR5(f,z) ) 

and FcR2(y:chan) = 

( new f@1.0:chan 

!fil(f)*2.0; FcR4(f,y) ) 

and FcR3(z:chan) = 

( new f@1.0:chan 

!fil(f)*2.0; FcR5(f,z) ) 



and FcR4(f : chan, y : chan) = 

( ?phosphz3(z) ; FcR7(f,y,z) ) 



let PhosphOO = 

( new phosph@l . : chan 

do !phosphy2(phosph)*l .0; 

Phosphl (phosph) 
or ! phosphz3 (phosph) *1 . 0; 
Phosphl (phosph) ) 

and Phosphl (phosph: chan) = 




run 1000 of FcROO 
run 1000 of IgGOO 
run 1000 of PhosphOO 



6 Appendix B 

site f on FcR associates site i on IgG with rate 2.0 

site y on FcR gets phosphorylated if site s on FcR is bound 

site z on FcR gets phosphorylated if site s on FcR is bound 

site s on FcR associates site sr on Src if site f on FcR is bound 

site s on FcR dissociates site sr on Src 

The SPiM code resulting from the automated translation of this model. 



directive sample 10.0 
directive plot FcR15(); 



FcR14(); FcR13(); FcR12(); 
FcRllO; FcRlOO; 
FcR9(); FcR8(); FcR7() ; 
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FcR6(); FcR5(); 
FcR4(); FcR3(); FcR2(); 
FcRK); FcROO; 
IgGK); IgGOO; 
PhosphK); PhosphOO; 
SrclO; SrcOO 

new f il@l . : chan(chan) 
new phosphx2<31 . : chan(chan) 
new phosphy3@l . : chan(chan) 
new ssr4@l . : chan(chan) 
new nil@0 . : chan 



let FcROO = 
( new f@1.0:chan 

!fil(f)*2.0; FcRl(f) ) 

and FcRl(f :chan) = 
( new s 100. 50: chan 

!ssr4(sl)*1.0; FcR5(f,sl) ) 

and FcR2(sl:chan) = 
( new f@1.0:chan 

do !fil(f)*2.0; FcR5(f,sl) 
or ?phosphx2(x) ; FcR8(sl,x) 
or ?phosphy3(y) ; FcR9(sl,y) 
or !sl; FcROO or ?sl; FcROO ) 

and FcR3(x : chan) = 
( new f@1.0:chan 

!fil(f)*2.0; FcR6(f,x) ) 

and FcR4(y:chan) = 
( new f@1.0:chan 

!fil(f)*2.0; FcR7(f,y) ) 

and FcR5(f : chan, si : chan) = 
( do ?phosphx2(x) ; FcRlKf ,sl,x) 
or ?phosphy3(y) ; FcR12 (f , si ,y) 
or !sl; FcRl(f) or ?sl; FcRl(f) ) 

and FcR6(f : chan, x : chan) = 
( new s 1<§0. 50: chan 

!ssr4(sl)*1.0; FcRll (f , si ,x) ) 

and FcR7(f : chan, y : chan) = 
( new s 1<§0 . 50 : chan 

!ssr4(sl)*1.0; FcR12(f , si ,y) ) 

and FcR8 (si : chan, x: chan) = 
( new f@1.0:chan 

do !fil(f)*2.0; FcRll (f , si ,x) 
or ?phosphy3(y) ; FcR14(sl ,x,y) 
or !sl; FcR3(x) or ?sl; FcR3(x) ) 

and FcR9 (si : chan, y: chan) = 



( new f (§1.0: chan 

do !fil(f)*2.0; FcR12 (f , si ,y) 
or ?phosphx2(x) ; FcR14(sl ,x,y) 
or !sl; FcR4(y) or ?sl; FcR4(y) ) 

and FcRlO (x : chan, y: chan) = 
( new f (§1.0: chan 

!fil(f)*2.0; FcR13(f,x,y) ) 

and FcRll (f : chan, si : chan, x : chan) = 
( do ?phosphy3(y) ; FcR15(f , si ,x,y) 

or !sl; FcR6(f,x) or ?sl; FcR6(f,x) ) 

and FcR12 (f : chan, si : chan, y: chan) = 
( do ?phosphx2(x) ; FcR15(f , si ,x,y) 

or !sl; FcR7(f,y) or ?sl; FcR7(f,y) ) 

and FcR13 (f : chan, x : chan, y: chan) = 
( new sl@0.50:chan 

!ssr4(sl)*1.0; FcR15(f , si ,x,y) ) 

and FcR14 (si : chan, x : chan, y: chan) = 
( new f (§1.0: chan 

do !fil(f)*2.0; FcR15 (f , si ,x ,y) 

or !sl; FcR10(x,y) or ?sl; FcR10(x,y) ) 

and FcR15 (f : chan, si : chan, x : chan.y: chan) = 
( do !sl; FcR13(f,x,y) or ?sl; FcR13(f,x,y) ) 



let IgGOO = 

( ?fil(i); IgGl(i) ) 

and IgGl(i:chan) = 
O 



let PhosphOO = 

( new phosphOl . : chan 

do !phosphx2(phosph) *1 .0; 
Phosphl (phosph) 
or ! phosphy3 (phosph) *1 .0; 
Phosphl (phosph) ) 

and Phosphl (phosph: chan) = 
O 



let SrcOO = 

( ?ssr4(srl); Srcl(srl) ) 

and Srcl (srl : chan) = 

( do !srl; SrcOO or ?srl; SrcOO ) 



(* run 1000 of ... *) 



